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ABSTRACT 

We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the 
Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally 
developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find 
isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorpo- 
rates several advances, including a metric-based gridding of the search parameter space (frequency, frequency 
derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies 
between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs J 1803-2 149 
and J21 1 1+4606, are young and energetic Galactic-plane pulsars (spin-down power above 6x 10 35 erg s 1 and 
ages below 100 kyr). The seven remaining pulsars, PSRs J01 06+4855, J0622+3749, J 1620-4927, J1 746-3239, 
J2028+3332, J2030+4415, J2 139+47 16, are older and less energetic; two of them are located at higher Galac- 
tic latitudes (\b\ > 10°). PSR J0106+4855 has the largest characteristic age (3 Myr) and the smallest surface 
magnetic field (2x 10 n G) of all LAT blind-search pulsars. PSR J21 39+47 1 6 has the lowest spin-down power 
(3x 10 33 erg s~‘) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency 
observations, only PSR J0106+4855 has detectable pulsations in the radio band. The other eight pulsars belong 
to the increasing population of radio-quiet gamma-ray pulsars. 
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1. INTRODUCTION 

The Fermi Gamma-ray Space Telescope has been operat- 
ing since launch in 2008 June. The Large Area Telescope 
(LAT) on board the Fermi satellite has an effective area of 
~ 0.8 m 2 (on-axis, above 1 GeV) and is sensitive to photons 
with energies from 20 MeV to more than 300 GeV (Atwood 
et al. 2009). Post-processing assigns an arrival time, energy E, 
and direction to the photons; we call these “events”. The ar- 
rival times are accurate to better than 1 ps, and the energies 
are accurate to within 15% between 0.1 and 10 GeV on-axis. 
The directional precision is energy-dependent: 68% of pho- 
tons have angular offset less than ~ 0?8 x(E/ GeVp*' 8 from 
the true direction (Abdo et al. 2009c). 

Gamma-ray pulsars are among the most interesting sources 
observed by the Fermi LAT. These are rapidly-spinning neu- 
tron stars whose regular beam of gamma-ray emissions passes 
by the detector with each rotation. The Fermi LAT has de- 
tected gamma-ray pulsations from more than 50 normal and 
millisecond pulsars (MSPs) previously discovered in radio- 
frequency searches (see e.g.; Abdo et al. 2010c; Ransom et al. 
2011; Keith et al. 2011; Coanard et al. 201 1; Theureau et al. 
2011 ). 

In contrast, so-called “blind” searches for gamma-ray pul- 
sars are not guided by any prior knowledge (e.g. from ra- 
dio or X-ray observations) of the pulsars’ parameters. Pre- 
vious blind searches of the data recorded by the Fermi - LAT 
have been spectacularly successful; they have discovered 26 
gamma-ray pulsars (Abdo et al. 2009a; Saz Parkinson et al. 
2010, 2011). This paper describes the “blind” discovery of 9 
additional pulsars using a new method. 

The blind search problem is computationally demand- 
ing because the relevant pulsar parameters (typically sky- 
position, frequency /, and spin-down rate /) are not known a 
priori and must be explicitly searched (Chandler et al. 2001). 
So far, blind searches of LAT data have used a clever “time- 
differencing technique” as described in Atwood et al. (2006); 
Ziegler et al. (2008). One powerful motivation for seeking 
even better methods is the application to MSPs, of which pre- 
vious blind searches have not found any. The blind search for 
MSPs is more computationally challenging due to the higher 
frequency range that must be covered. Moreover, most MSPs 
are in binary systems, where the orbital modulation parame- 
ters must also be searched, adding orders-of-magnitude to the 
computational complexity and challenge. 

This paper presents results from a new effort to find isolated 
gamma-ray pulsars (including MSPs), using a novel method 
inspired by computationally-efficient techniques recently de- 
veloped to search gravitational-wave detector data for weak 
continuous-wave signals from rapidly-spinning isolated neu- 
tron stars (Pletsch & Allen 2009; Pletsch 2010). In partic- 
ular, the search method uses a recently-developed optimal 
incoherent-combination method first described in Pletsch & 
Allen (2009) together with a "sliding coherence window tech- 
nique” (Pletsch 2011). 

The method was originally intended to find isolated MSPs 
up to 1.4 kHz spin frequency with spin-down rates in the 
range -5x 10~ B Hz s" 1 < / < 0, with a corresponding range 
of characteristic ages r = -f/2f. However, the search tech- 
nique is also sensitive to normal (non-MSP) isolated pulsars, 
and several normal pulsars were discovered soon after we be- 
gan. This prompted us to extend the spin-down-rate search 
range, in order to include younger objects (having smaller r), 
down to ages (~ kyr) comparable to the Crab pulsar. Here 


we report on the discovery of nine pulsars from this ongoing 
effort in which a total of 1 09 sources selected from the Ferrni- 
LAT Second Source Catalog (Abdo et al. 2011b) are being 
searched for new gamma-ray pulsars. 

The outline of this paper is as follows. Section 2 describes 
the LAT data preparation, and the selection of unidentified 
sources to search for previously unknown gamma-ray pulsars. 
The new hierarchical search method is explained in Section 3 
and illustrated in Section 4 with a detailed example: our first 
discovery, PSR 11620^1927. Section 5 presents the results for 
all the new gamma-ray pulsars. The search for counterparts in 
other regions of the electromagnetic spectrum is discussed in 
Section 6. This is followed by a brief conclusion. 

2. SOURCE SELECTION AND DATA PREPARATION 

The Fermi-LAI Second Source Catalog (2FGL, Abdo et al. 
2011b) lists 1873 sources, described by fits to elliptically- 
shaped 95%-confidence sky regions. Among these sources 
576 are not associated with counterparts observed at other 
wavelengths and thus might contain unknown gamma-ray 
pulsars. 

In searching for new gamma-ray pulsars, it is important 
to identify and exclude sources that are blazars. Previously- 
observed Fermi-hKT gamma-ray pulsars (see e.g. Abdo et al. 
2010c) have sharp cutoffs in their emission spectra at a few 
GeV, and stable gamma-ray fluxes. In contrast, blazars emit 
above 10 GeV, and their fluxes vary with time. An illustra- 
tion of the different spectral properties and variability behav- 
ior can be seen in Figure 17 of Abdo et al. (201 lb). Gamma- 
ray pulsars tend to have large curvature significances (“Sig- 
nif_Curve” in the 2FGL catalog, which gives the improve- 
ment in the quality of the spectral fit when changing front a 
power law' to a curved spectral model) and small variability 
indices (“Variability _Index” in the 2FGL catalog, a measure 
of flux instability over time). In contrast, blazars tend to have 
small curvature significances and large variability indices. 

For the search for gamma-ray pulsations we select 2FGL 
sources with curvature-significance values greater than 4a 
(here and throughout the manuscript, a denotes the stan- 
dard deviation for a Gaussian distribution), and variability in- 
dices smaller than 41.6. We further select bright objects by 
choosing sources with detection significances (“Signif_Avg”) 
greater than 10a. Finally, we restrict the list to sources with 
no known associations (“unassociated sources”) or associated 
with known supernova remnants (SNRs). By applying the 
above selection criteria to a preliminary version of the 2FGL 
catalog, we obtained a list of 109 2FGL sources to search 
for gamma-ray pulsations. (One source, 2FGL J062 1 .9+3750 
was associated with an active galactic nucleus (AGN) in the 
final 2FGL catalog.) 

The Fermi Science Tools (ST) 25 v9r23pl are employed to 
select the Fermi - LAT events for our search. Using gtselect , 
we take events from 2008 August 4 to 201 1 April 6, with re- 
constructed directions within 8° of the gamma-ray sources, 
energies above 100 MeV, and zenith angles < 100°. Only 
events belonging to the Pass 6 "Diffuse” class are retained, 
as those events have the highest probability of being photons 
(Atwood et al. 2009). For the event probability weighting de- 
scribed below we use the P6_V 1 1 Instalment Response Func- 

25 See http://fermi.gsfe.nasa.gov/ssc/data/ 
analysis/scitools/overview . htral and http: //fermi . 

gsf c . nasa . gov/ ssc/data/analysis/documentation/ 
Cicerone/ for details and tutorials. 
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tions (IRFs). Using the gtmktime tool, times when the rock- 
ing angle of the satellite exceeded 52° are excluded. We also 
require that DATA_QUAL and LAT_CONFIG are unity, and 
that the Earth’s limb does not impinge upon the Region Of 
Interest (ROI). 

As demonstrated in previous work (Bickel et al. 2008; Kerr 
201 la), the sensitivity of gamma-ray pulsation searches can 
be improved by weighting the jth event with a probability 
wj e [0, 1] that it originated from a putative pulsar. The pho- 
ton weights Wj are calculated by using a full spectral model 
of the region around the gamma-ray source, and by exploiting 
the IRFs to provide background rejection which is superior to 
simple angular and energy data-selection cuts. For the first 
time in this paper such a photon-weighting scheme has been 
applied in a blind search; further details are given in Section 3. 

To calculate the weights wj for each event from a given 
source, we perform likelihood spectral analyses using the 
pointlike tool (see Kerr 2011b, for a description). For each se- 
lected source, we construct a spectral model for the region by 
including all sources of the 2FGL catalog found within 8° of 
the selected source, using the spectral forms given in the cat- 
alog. The spectra of the selected sources are modeled as ex- 
ponentially cut-off power laws, typical of known gamma-ray 

pulsars, of the form No (E / GeV) 1 exp (~E /E c ) , where N 0 is 
a normalization factor, T is the photon index and E c is the cut- 
off energy. The Galactic and extragalactic diffuse emission 
and residual instrument background also enter the calcula- 
tion of the weights. The Galactic diffuse emission is modeled 
using the gllJem_v02_P6_Vl 1 _D1FFUSE map cube, while 
the extragalactic diffuse and residual instrument backgrounds 
is modeled using the isotropic _iem_v02_P6_V 11 JDIFFUSE 
template (a detailed description of these background models 
can be found in Section 3 of Abdo et al. (2010a)). These 
models are available for download at the Fermi Science Sup- 
port Center 26 . The tool gtsrcprob is then used to calculate 
the event weights wj based on the best-fit spectral models ob- 
tained from the maximum likelihood analyses. 

3. THE NEW SEARCH METHOD 

In a year, the LAT detects of order 1 0 3 photons from a typ- 
ical gamma-ray pulsar; in the same year, a typical pulsar ro- 
tates at least 10 s times around its axis. The blind-search prob- 
lem is to find a rotational-phase model <I>(f) = 2rr(/r -t-/V 2 /2) 
and a sky-position that match the Solar System Barycen- 
ter (SSB) arrival times t of the different photons, where <I> 
denotes the rotation angle of the star about its axis, in ra- 
dians, measured from its starting position at t = 0, and ob- 
served at the SSB. The signal hypothesis is that the photons’ 
arrival times are “clustered” near specific “orientations” of the 
star (i.e., $(/) mod 2~ deviates from uniformity on the inter- 
val [0,27r]). The null hypothesis is that the arrival times of 
the photons are a random Poisson process. In this paper we 
do not explicitly indicate the dependence of $ on f, / and 
sky position, but this dependence is important and implicit in 
many formulae below. 

To find a matching phase-model, a grid of “templates” in 
the four-dimensional parameter space of sky position and 
(/,/) is constructed. Note that the 2FGL catalog sky posi- 
tions of the targeted unassociated sources based on the spatial 
distribution of events are typically not precise enough for pul- 
sar searches. A search grid of sky points around this catalog 

26 http://fermi.gsfc.nasa.gov/ssc/data/access/Iat/BackgroundModels.html 


position is needed to reduce signal loss arising from imperfect 
correction of the Doppler shifts caused by the Earth’s orbital 
motion around the SSB. The need for sky gridding is partic- 
ularly acute for MSP spin frequencies. Therefore, in contrast 
to previously published blind searches 27 , we grid a circular 
sky-region centered on the 2FGL catalog source location us- 
ing a radius which is 20% larger than the semi-major axis 
of the 95% confidence elliptical error region (given by the 
“Conf_95_SemiMajor” parameter). 

Unfortunately the number of templates (grid points) re- 
quired to discretely cover the entire four-dimensional search 
parameter space increases as a high power of the coherent 
integration-time. Hence a fully-coherent approach for several 
years of data is computationally impossible. Therefore, we 
employ a search strategy which is designed to achieve maxi- 
mum overall sensitivity at fixed computing cost 28 . 

To efficiently scan through years of Fermi-LAT data for 
previously unknown gamma-ray pulsars, we use a so-called 
hierarchical search approach. This is analogous to hierarchi- 
cal methods used in searches for gravitational-wave pulsars 
(Schutz & Papa 1999; Papa et al. 2000; Brady & Creighton 
2000; Abbott et al. 2009a,b; Cutler et al. 2005). In a first 
semi-coherent stage, we here adopt the optimal metric-based 
gridding methods described in Pletsch & Allen (2009) along 
with the sliding coherence window technique (Pletsch 2011). 
In a second stage, significant semi-coherent candidates are au- 
tomatically followed up in a fully-coherent analysis. Finally, a 
third stage further refines coherent pulsar candidates by using 
higher harmonics. Full details of the complete search scheme 
will be presented in forthcoming work (Pletsch & Guillemot 
2011 ). 

Here we first describe the principle of the method, to firmly 
establish the analogy with the existing gravitational-wave lit- 
erature. Then we describe what is done in practice, which is 
mathematically equivalent (up to justifiable approximations) 
but computationally more efficient. 

In the first stage, a semi-coherent detection statistic S is 
computed for each template. We refer to S as “semi-coherent” 
because it is effectively the incoherent sum over several years, 
of terms which are coherent over several days. The coher- 
ent terms are the power in Fourier bins, calculated by treating 
each photon arrival as a delta function in time. 

Denoting the arrival time of the /th event (photon) at the 
SSB by tj, the coherent power P r in a (Gaussian) window 
centered at time r is defined by 


Pr = 


;=i 


-<'4>(Z j) j-r) 2 /T 2 


(1) 


The sum is taken over all photons (here, N = 8000) in the data 
set; the effective window duration is f e~ 2lTT " :T 'dr = T j\fl. 
As described in Section 2, the weights Wj estimate the proba- 
bility that the photon comes from the selected source. 

To form the semi-coherent detection statistic S. the values 


27 Ziegler et ai. (2008) argues correctly that for two-week data-stretches, 
sky gridding is not essential. However for data stretches of length compara- 
ble to a year or longer, sky gridding is necessary to avoid significant loss of 
signal-to-noise ratio. 

28 If this constraint is removed, then it is obvious that a more sensitive 
method exists. 
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of P r are summed (“incoherently combined”) 



Note that in this definition we have subtracted a constant 
(phase model indepdendent) term. Because it contains a 
Gaussian window, the integrand in Equation (2) falls off ex- 
ponentially at early and late time (large values of |r|). Thus 
the limits of integration can be taken as the entire real line 
r € (—oo, oo); to good approximation this gives the same value 
as integrating only over the total observation interval (about 
975 days in this search). 

The semi-coherent detection statistic S is an incoherent sum 
of powers, which discards the phase information overtime pe- 
riods longer than of order T. This uniform overlap maximizes 
the search sensitivity for fixed T and for fixed computational 
resources (Pletsch 2011). For computational efficiency, in this 
search we choose the N = 8000 photons with the highest prob- 
abilities (largest values of wj). 

To understand how to compute S efficiently, one can ex- 
plicitly evaluate Equation (2). Completing the squares in the 
product of the Gaussians and carrying out the integration over 
r, one obtains 

S = ]C WjWke~ imj) ^ m e~* (,J - tl)2/T2 - ]T w 2 . (3) 

7=1 tel 7=1 

Here the effective duration of the Gaussian window is 
f e~' KT ' /r dr = T. In practice, to compute the semi-coherent 
power efficiently, we replace the Gaussian window in Equa- 
tion (3) with a rectangular window of the same duration 7\ as 
given below in Equation (5). In this search, the width of the 
rectangular window is T = 2 19 s (« 6 days). 

The template grid in parameter space is the Cartesian prod- 
uct of a rectangular two-dimensional grid in (/,/) and a sky- 
grid which has constant density when orthogonally projected 
onto the ecliptic plane. The problem of constructing efficient 
search grids has been intensively studied in the context of 
gravitational wave searches (see e.g. Brady et al. 1998; Brady 
& Creighton 2000; Prix 2007; Pletsch & Allen 2009; Pletsch 
2010) and we employ these concepts here. The values of fre- 
quency are equally spaced, separated by the FFT frequency- 
bin width Af = 1 /T. In practice, to reduce the fractional loss 
in S for frequencies not coinciding with Fourier frequencies, 
we use a computationally-efficient interpolation, referred to 
as “interbinning” (Ransom et al. 2002). The spacing in the 
other three dimensions is determined by a metric which mea- 
sures the fractional loss in the expected value of 5, that arises 
if the signal is not located exactly at a grid point (Balasubra- 
manian et al. 1996; Owen 1996; Prix 2007). In spin-down /, 
we use a uniform grid spacing Af = s/llOm/tt^T 2 . In this 
search m = 0.3 is the maximum tolerable fractional loss in S, 
and from Pletsch & Allen (2009), 



where t = Jf-f/N is the mean photon arrival time. The de- 
scription of this grid can be found in Pletsch & Allen (2009) 
with a detailed derivation in Pletsch (2010). The grid in the 
sky is determined by the same metric, permitting a maxi- 


mum fractional loss in in the value of S. The spacing of the 
sky grid is determined by the Doppler-shift arising from the 
Earth’s (more precisely, the Fermi satellite’s) motion around 
the Sun. At the north Ecliptic pole the angular spacing is 
Ad = Vlmc/vfD , where c is the speed of light and D is a 
baseline distance (defined below). When the entire sky grid 
is projected into the plane of the ecliptic, the grid-points are 
uniformly spaced on the plane (Astone et al. 2002; Abbott 
et al. 2009a,b). This angular spacing is similar to the an- 
gular spacing in the diffraction pattern of a two-slit system, 
where the wavelength is c/f and the separation of the two 
slits is the straight-line distance D between two points on the 
Earth's orbit about the Sun. If the coherent integration time T 
is less than half a year, then D = (998 s ) c sin(~T / 1 yr). If the 
coherent integration time T is greater than half a year, then 
D = (998 s)c is the diameter of the Earth’s orbit about the Sun. 

To compute S efficiently, a time series is constructed and 
subsequently Fourier-transformed into the frequency domain. 
The time series contains T A / BW bins, where A / BW is the to- 
tal frequency bandwidth being searched at a time using com- 
plex heterodyning 29 at the center of A / BW (see e.g. Patel et al. 
2010). The time series is initialized to zero, then the values 
of are added into the bins determined by the 

time differences A t# - tj-t a, for all pairs of photons j,k for 
which 0 < \Atjk\ < T ; the bin index is obtained by round- 
ing the absolute value of the product A / BW Atjk to its nearest 
integer value. 

Then the array is Fourier-transformed into the frequency 
domain (exploiting the FFT) to obtain S over the entire f grid. 
Up to an overall window-dependent normalization, one can 
write for given values of /, / and sky position, 

N 

S = ]T G(| Af;*|/r ) wj w k e- 2 ” fA (5) 

7, tel 

where the rectangular function Q(x ) is unity if 0 < x < 1 and 
vanishes otherwise 30 . 

Although other aspects are different, the use of an FFT ap- 
plied to time differences is very similar to techniques pre- 
viously used in blind searches of Fermi -LAY data (Atwood 
et al. 2006, Equation 3 therein has a typo which is corrected 
in Equation 2 of Ziegler et al. (2008)). The method used in 
Atwood et al. (2006) was the first application of this classic 
method (e.g. Blackman & Tukey 1958) to gamma-ray astron- 
omy (in estimating the power spectrum an approximate auto- 
covariance function is calculated using a maximum lag, and 
then Fourier-transformed). 

In contrast to previous searches, our method uses an optimal 
gridding of the parameter space for both the semi-coherent 
and coherent stages, as well as an automated follow-up, and 
incorporates the spin-down corrections in a way that permits 
heterodyning and highly-efficient code. 

The search was done on the 1680-node Atlas Computing 
Cluster (Aulbert & Fehrmann 2008) built around four-core 
processors with 8 GB of random-access memory; for these 
we used a heterodyning bandwidth A / BW = 256 Hz. Breaking 
the full frequency- range of the search into frequency-bands 
allows the computation to fit into memory, and also allows the 

29 Complex heterodyning is a procedure which shifts frequencies in time- 
series data by a fixed offset /),. This is accomplished by multiplying the time 
series by e~ 2 w f hi J, shifting all frequencies by f h . 

30 Note that by symmetry S is real, because interchanging indices j and k 
is equivalent to complex conjugating the exponential factor. 
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Fig. 1 . — Search results for the newly-discovered pulsar PSR J 1 620-4927 (the large black dot in the lower-left of each panel). The left panel (a) shows the 
semi-coherent search results, representing about 2 CPU-years of computing on a single core. (Since the computing cost scales as the square of frequency, a search 
up to 64 Hz as in previous searches (Abdoet al. 2009a: Saz Parkinson et al. 2010) would have taken about 1.4 CPU-days on a single core.) The bottom-left panel 
shows the semi-coherent detection statistic S as a function of / and /, maximized over the sky-grid. The value of S is represented by the colorbar. A further 
maximization over f is shown to the right; a further maximization over / is shown above. PSR J 1 620-4927. the darkest point near the bottom left of each figure, 
stands out clearly from the noise. In the same form, the right panel (b) presents the fully-coherent follow-up search results of the previous candidate (and every 
other candidate “dot”) shown in the left panel. The quantity plotted is now the fully-coherent detection statistic P over the entire data set. As explained in the 
text, for each candidate this covers a region of parameter space which is four steps of the semi-coherent grid in each dimension. 


use of different sky grids in each band. This further reduced 
the computational cost, since the number of required sky grid 
points increases with the square of frequency. 

After computing S on the four-dimensional grid in param- 
eter space, points with statistically-significant values of S are 
candidates for possible pulsar signals, and are followed up in 
a second stage. This is done by “refining the grid” and in- 
creasing the coherent integration time. This is an hierarchical 
scheme which is analogous to “zooming”: successively swap- 
ping microscope objectives for ones of higher magnification, 
then re-centering the interesting point on the slide (see e.g. 
Cutler et al. 2005; Krishnan et al. 2004). In our case, this is 
done by constructing the fully-coherent detection statistic P 
over the entire data set (or equivalently taking the Gaussian 
window-size T —> oc in Equation (1)) obtaining 


P = 


E w -> e 

7=1 


■*&(<>) 


( 6 ) 


where for convenience we have normalized by the positive 
constant k given by 


E 

7=1 


ir: 


(7) 


The computing cost to coherently follow up a single candidate 
is negligible in comparison to the cost of the previous semi- 
coherent search. 

In selecting statistically-significant semi-coherent candi- 
dates which are automatically followed up in the second stage 
using a fully-coherent analysis, we do not use a fixed thresh- 
old to define “statistical significance”. In the semi-coherent 
stage, the search code keeps an internal list of the strongest 
signal candidates. Each member of this list is coherently fol- 
lowed up and corresponds to the largest value of S detected in 


eight adjacent spin-down values for the entire heterodyning 
frequency bandwidth and a single sky point. 

The refined grid of the fully-coherent follow-up covers a 
region of parameter space of size (4A /) x (4 A/) x (4A 0 x 
4A0) when projected into the ecliptic plane. In other words, 
it covers a region whose volume is 256 times larger than the 
volume of a fundamental cell in the original grid: its extent 
in each dimension of parameter space is four grid-intervals. 
The refined grid has a spacing given by the previous formulae 
for A/, A /, and Ad, except that the coherence time T is set 
equal to the length of the entire data set, and 7=1. Since in 
this case only a small parameter-space region around the can- 
didate is explored, it is computationally-efficient to compute 
P directly in the time-domain (FFTs are not used), exploiting 
the sparsity of the photon data. 

If the value of P, which measures the fully-coherent power 
(in a single harmonic), is statistically significant, then in a 
third stage further refinement is carried out using higher har- 
monics (Fourier-components). We adopt the so-called H- test, 
which has been widely used in X-ray and gamma-ray pulsar 
detection (de Jager et al. 1989; de Jager & Biisching 2010). 
This test measures the statistical significance of the energy in 
the first 20 (non-DC) Fourier-components of the pulse-profile 
as a function of phase. As in Equation (1), the H- test can 
also be modified to include the photon probability weights Wj 
(see Kerr 20 1 la). Note that Equation (5) in Kerr (201 la) con- 
tains an error; corrected formulae used in this work are given 
below. 

The weighted H - test statistic is defined as follows. For each 
photon arrival time tj, the pulse phase Xj (between zero and 
one) is calculated as Xj = | <I>(t ; ) mod 2rc] /2 tt. The pulse pro- 
file (for 0 < .v < 1 ) is a sum over the photons 

N 

p(x) = WjS(x-Xj), 


(8) 
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Fig. 2. — The weighted H - test statistic for PSR J 1620-4927. Panels (a) and (b) show contour plots of the weighted H - test statistic as a function of sky 
position (a) and (/,/) (b). These are peaked at the parameter-space location indicated by the black cross and the axes show the offset from these values. Panel (c) 
shows how the maximum weighted H - test statistic [shown in (a) and (b)] accumulates with time. The //-test increases (approximately) linearly with time, as 
expected for a pulsar that is emitting uniformly. 


where Xj is the pulse phase of the / th photon and d(x) is a 
one-dimensional Dirac delta-function. It can be expressed as 
a Fourier series 

OG 

/?(x) = k a«e 27ri£ \ (9) 

£~~oo 

which implies that the (complex) Fourier coefficients are 
given by 

1 

Ct£ — — 

K 

where the definition of k is identical to Equation (7). The nor- 
malization of the Fourier coefficients has been chosen so that 
if the photon arrival times are uniformly distributed, indepen- 
dent random variables, then in the limit of large numbers of 
photons, for £ > 0, R(oy) and Q ; (o/>) are independent Gaussian 
random variables with zero mean and unit variance. Finally, 
the weighted H- test statistic is defined as 


z 


Wit 


-2-7rtf'Xj 


( 10 ) 


H = 


max 

I <£,<20 



— 4(L- 1) . 


(11) 


The quantity which is subtracted, 4(L- 1), is motivated by an 
empirical numerical study (de Jager et al. 1 989), providing the 
best omnibus test for unknown pulse profiles. 

Maximizing H over sky-position, frequency /, and spin- 
down rate / selects the “narrowest and sharpest” overall pulse 
profile. In contrast, maximizing the statistic / > =|qi| 2 of 
Equation (6) favors putting more power into the lowest har- 
monics. 

Using H as the test statistic, a further stage of parameter- 
space refinement is done in the same way as before: in the fre- 
quency, spin-down, and sky parameters, we cover regions that 
include four grid-steps (in each dimension) of the previous 
grid. The chosen refinement factor (ratio of the number of grid 
points after and before refinement, in each dimension) is about 
an order of magnitude. At each grid point, the weighted H- 
statistic of Equation (11) is found. The parameter-space point 
with the largest statistic is selected as our best estimate of the 
pulsars’ parameters which are then further refined through the 


timing-analysis procedure (Ray et al. 201 1) described in Sec- 
tion 5.2. 

The hierarchical search pipeline has been validated by suc- 
cessfully recovering previously known gamma-ray pulsars, 
including some of the brightest gamma-ray MSPs (Abdo et al. 
2010c; Pletsch & Guillemot 2011). In the next section, the 
complete search scheme is illustrated with a detailed exam- 
ple. 


4. EXAMPLE: RESULTS FOR PSR J 1620-4927 

Before giving the results for all of the pulsars that have 
been discovered with this new search method, we first go 
through a single example in detail. This illustration uses 
PSR J 1620-4927, the first new pulsar found in this work. 

Figure 1 (a) shows the first stage of the analysis. For each 
/ and f value, the largest value of S found in the sky grid 
is displayed as an intensity. The point of highest intensity 
corresponds to PSR J 1 620-4927’s f and / parameters. 

The second analysis stage is the automated follow-up of all 
candidates shown in Figure 1(a). As explained in Section 3, 
this is accomplished by carrying out a fully-coherent search 
over a small region of parameter space around each candi- 
date which has been identified as statistically significant in 
the previous semi-coherent stage. If the candidate found in 
the semi-coherent step was simply a statistical outlier, then 
the fully-coherent statistic P will not be significant. Fig- 
ure 1(b) presents the results of the fully-coherent statistic P for 
PSR J 1 620-4927. Again the point of highest intensity com- 
pared to the background, indicating the presence of a coherent 
signal in the data set, is due to the new pulsar. 

The importance of the photon probability weights wj can 
be illustrated by repeating the analysis using the same 8000 
photons with the weights set to unity. (Note that the weights 
were used in selecting the 8000 photons, so this is not a 
complete comparison.) The result is that Figure 1(a) still 
shows a statistically-significant outlier in the semi-coherent 
search step, which is followed up automatically and gives a 
statistically-significant outlier in the fully-coherent output of 
Figure 1(b). In both steps the signal and its statistical sig- 
nificance are reduced, but the pulsar is still found. However 
for some of the other new pulsars, this is not the case: if the 
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weights are set to unity, then the pulsar is not detected. Over- 
all the weights play a larger role for sources in crowded re- 
gions of the sky, for example near the Galactic plane. 

At the third stage, further refinement is carried out by max- 
imizing the weighted H- test statistic over a small region of 
parameter space around the most significant candidate from 
the previous step. As described in Section 3, the parameter- 
space grid used at this stage is yet another order-of-magnitude 
finer than the one of the previous stage. Figure 2 shows the 
corresponding results for PSR J 1620-4927. 

The parameter-space point with the largest H- test statistic 
is selected as our best estimate of the pulsar parameters at this 
stage. The uncertainties in these estimated parameters can be 
obtained by using the fact that a lex deviation for a Gaussian 
distribution has a value m 0.6 of the maximum. In the current 
example, as shown in Figure 2, the value of the maximum 
weighted /7-test statistic is 566. 

Figure 2(c) illustrates the way that the weighted H- test 
statistic accumulates over the total observation time. The lin- 
earity of this plot indicates that the properties of the source, 
and of the data, were time-stationary; it shows that the pul- 
sar’s emission is not changing with time. In this way it also 
gives additional confidence that the pulsar is real and not a sta- 
tistical noise outlier. (We have not provided analogous plots 
for the other eight pulsars reported in this paper, but they are 
similar.) 

As shown in de Jager & Biisching (2010) and Kerr (201 1 a), 
the weighted H - test has a false alarm probability approxi- 
mately described by P? A ~ e A)AH . For this example, the as- 
sociated false alarm probability is log K) fp A « -98. We have 
not tried to rigorously estimate the trials factor, but since 
the total number of floating-point computations that could be 
performed by our computing systems over a period of some 
months is less than 10 21 , we conclude that PSR J 1 620— 4927 
is a real gamma-ray pulsar and not a statistical outlier. 

5. RESULTS FOR ALL NEW PULSARS 

Table 1 shows the names and sky positions of the nine 
discovered pulsars, and also lists known source associations. 
The inferred rotational parameters of the new pulsars are pre- 
sented in Table 2. In addition, further derived parameters are 
given, including the spin-down power E = -At r 2 / //, where 
I = 10 45 g cm 2 is the assumed neutron star moment of iner- 
tia, and the magnetic field strength at the neutron star sur- 
face B s ~ 3.2 x 10 ,9 (-/// 3 s- 1 ) 1 ' 2 G and at the light cylinder 
Bi c ~ 2.94 x 10 8 (-// 3 s 5 ) 1,2 G, respectively. 

Figure 3 plots the newly discovered pulsars on a f-f di- 
agram, where they can be compared with the known pulsar 
population. One can see that the new pulsars lie in a simi- 
lar region (and hence belong at large to the same population) 
as the previously-found blind-search Fermi- LAT gamma-ray 
pulsars. 

Two of the nine new' systems, PSRs J 1803-2 149 
and J2 1 1 1 +4606, are young energetic pulsars 

(E > 6 x 10 35 erg s -1 and r<100kyr), located near the 
Galactic plane. Among the seven remaining less-energetic 
and older pulsars, five are also located near the Galactic plane 
(| b\ < 5°) and two are found at higher Galactic latitudes 
(|£>| > 10°). One of these, PSR JO 106+4855, has the largest 
characteristic age r = 3 Myr and the lowest surface magnetic 
field strength (B s « 2 x 10 n G) of all LAT blind-search 
pulsars found to date. Also standing out, PSR .12139+4716 
has the smallest spin-down power (£ = 3 x 10 3 , erg s~ ! ) 



FIG. 3. — The frequency / and spin-down rate / of the pulsar population. 
The black dots show the approximately 1800 pulsars in the ATNF catalog 
(Manchester et al. 2005). excluding pulsars in globular clusters. The green 
dots show the 26 gamma-ray pulsars discovered in previous blind Fermi - LAT 
searches (Abdo et al. 2009a; Saz Parkinson et al. 2010. 201 1). The nine red 
stars show the newly-discovered gamma-ray pulsars reported in this paper. 
The dotted lines indicate constant characteristic ages r, the dashed lines show 
contours of constant spin-down power E and the dashed-dotted lines signify 
contours of constant surface magnetic field strength B s . 


among all non-recycled gamma-ray pulsars ever detected (cf. 
De Luca et al. 2011). 

5.1. Source associations 

As listed in Table 1, all but one of the new pulsars have 
gamma-ray counterparts from the Fermi -LAT First Source 
Catalog (1FGL, Abdo et al. 2010a). The exception is 
PSR J2028+3332, which first appeared in the 2FGL catalog. 
In addition, these eight 1FGL unassociated sources have been 
classified as likely pulsar candidates based strictly on a com- 
parison of their gamma-ray properties with previously de- 
tected LAT pulsars (Abdo et al. 201 la). 

PSR J21 11+4606 is associated with the source 
0FGL J21 10.8+4606 from the Fermi Bright Source List 
(Abdo et al. 2009b). Only ten objects from this catalog still 
remain to be identified. 

The new system PSR J 1620-4927 is associated with the 
unidentified AGILE source 1 AGL J 1 624-4946 (Pittori et al. 
2009). This AGILE source is however about ten times 
brighter, with a flux of (67 ± 13) x 10~ 8 ph cm' 2 s~ ! . In ad- 
dition, it has an error radius of 0358 and encloses other 2FGL 
sources. The new gamma-ray pulsar could therefore only 
explain a fraction of the total flux detected by AGILE from 
1 AGL .11624-4946. 

While not formally associated in the 2FGL catalog (Abdo 
et al. 2011b). PSR J2028+3332 lies within the 99% error 
contours for the EGRET source 3EG J2027+3429. How- 
ever, the pulsar accounts for only about a fifth of the 
flux for the EGRET source. Another nearby LAT source, 
2FGL J2025. 1+3341, is associated with an AGN and has a 
peak LAT flux that can account for the remainder of the mea- 
sured EGRET flux. 

We have searched the TeVCat online catalog of TeV 
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sources 31 for very high-energy counterparts to the newly- 
discovered pulsars, but found no associations. This is not 
surprising, considering that most of the new gamma-ray pul- 
sars have moderate spin-down power (~ 10 34 erg s' 1 ) and rel- 
atively large characteristic ages (> 100 kyr), while pulsars 
associated with TeV Pulsar Wind Nebulae (PWNe) tend to 
be young and energetic. PSRs J 1803-21 49 and J21 1 1+4606 
have properties that make associations with TeV PWNe plau- 
sible; however, no detections of TeV emission from these di- 
rections have been reported. 

5.2. Timing analysis 

For each new gamma-ray pulsar the definitive parameters as 
given in Tables 1 and 2 are determined via precise measure- 
ments of pulse Times Of Arrival (TOAs) and by fitting the 
parameters of a timing model to these measurements using 
methods described in Ray et al. (201 1). 

To obtain a set of TOAs, photons are extracted for each 
source using a radius and minimum energy cut to optimize 
the signal-to-noise ratio for each pulsar. Then these data are 
subdivided into segments of about equal length. The set of 
best-guess pulsar parameters from the search (maximizing H) 
is used as a preliminary ephemeris to fold the photon ar- 
rival times, producing a set of pulse profiles. The TOAs are 
then measured by cross-correlating each pulse profile with a 
multi-Gaussian or kernel-density template derived from fit- 
ting the entire data set. This is done using the unbinned max- 
imum likelihood method described in Ray et al. (2011). Then 
Tempo2 (Hobbs et al. 2006) is used to fit the TOAs to a tim- 
ing model including sky position, frequency and frequency 
derivative. 

As often is the case, the youngest object found in this work 
(PSR J21 11+4606, r = 17 kyr), exhibits particularly large 
timing noise. In turn, this requires including higher-order 
frequency-derivative terms to make a good timing model fit; 
in this case up to the fourth derivative of the frequency as 
shown in Table 2. For the same reason, the timing models in 
Table 2 for PSRs J2030+4415 and J 1 803-2 1 49 include terms 
up to the third and second frequency derivative, respectively. 

Based on these timing solutions. Figures 4-12 show the 
resulting phase-time diagrams and pulse profiles for each of 
the newly-discovered pulsars. These plots are obtained from 
calculating the phase for each of the 8000 photons selected 
for the blind search and using the parameters listed in Ta- 
bles 1 and 2, and weighting each event by its probability of 
originating from the pulsar wj. The integrated pulse profiles 
(weighted pulse phase histograms) are constructed with a res- 
olution of 32 bins in phase per rotation. The 1 a error bars in 
the integrated pulse profiles are statistical only and are given 
by (EjW 2 j) 1/2 , where j runs over all events falling into the 
same phase bin. (Note that the formula for the fractional er- 
ror has the opposite sign of the exponential: -1/2.) 

As seen in Figures 4-12, eight of the nine pulsars have pulse 
profiles with two peaks. For pulsars where the two peaks 
are separated by nearly one-half of a rotation it is possible 
to detect or discover the pulsar at the second harmonic (i.e. 
at twice the actual spin frequency of the pulsar). For pulsar 
candidates which were discovered with (apparently) single- 
peaked profiles, we tested for the true fundamental spin fre- 
quency by folding at the subharmonic of the discovery fre- 
quency. If the subharmonic is the correct frequency, the two 

3 http : //tevcat . uchicago . edu/ 


resulting peaks may satisfy one or more of the following con- 
ditions: offsets that are measurably different than 0.5 in phase, 
significant differences in the integrated weighted counts un- 
der each peak, or significant differences in the shape of the 
peaks in different energy bands. With modest signal-to-noise 
ratios, such determinations are not always conclusive. For 
PSR J0106+4855, our identification of the true period was 
subsequently confirmed by the detection of radio pulsations 
with a single peaked profile. For the one apparently single- 
peaked pulsar in our sample, PSR J2 139+47 16, none of the 
above tests yield strong evidence for the profile being double 
peaked at half the frequency. Additional data will be required 
to strengthen this conclusion. 

In order to further characterize each pulse profile, we fit the 
pulsars’ weighted gamma-ray peaks to Lorentzian lines. The 
derived pulse shape parameters are listed in Table 3. Note that 
PSRs J0622+3749 and J 1746-3239 show indications of sub- 
structures in their main gamma-ray peak. For these pulsars 
a single Lorentzian line function is used for fitting the main 
component. Apart from PSR JO 106+4855, which is detected 
at radio wavelengths (cf. Figure 13), there is no particular 
reference for absolute phase of these pulsars. For the other 
eight pulsars we have arbitrarily assigned the absolute phase 
reference such that the first peak occurs at a value of 0. 1 in 
phase. The gamma-ray pulse profiles shown in Figures 4-12 
along with the pulse-profile parameters listed in Table 3 are 
very similar to those of the previously-discovered gamma-ray 
pulsars in blind searches (Abdo et al. 2009a; Saz Parkinson 
et al. 2010; Abdo et al. 2010c), further supporting the theory 
that the gamma-ray emission consists of fan beams produced 
in the outer magnetosphere. 


5.3. Spectral parameters 

The spectral parameters for the new pulsars are obtained 
by fitting each phase-averaged spectrum with an exponen- 
tially cut-off power-law with a photon index T and a cutoff 
energy E c . The results for each pulsar are listed in Table 4. In 
addition to T and E c which are explicit parameters of the fit, 
Table 4 also gives the important derived physical quantities of 
photon flux F \ oo (in units of photons cm -2 s' 1 ) and the energy- 
flux Gioo (in units of erg cm' 2 s' 1 ) for events with energies 
between 1 00 Me V and 1 00 GeV. 

Analogous to the pulse-profile properties, T and E c mea- 
sured for the new pulsars are also similar to those observed for 
previously-detected gamma-ray pulsars (Abdo et al. 2010c). 
This is not surprising, because (as described in Section 2) tar- 
get sources for our search have been selected based on similar- 
ity of their spectral properties to known gamma-ray pulsars. 

The distance (3 kpc) for one of the new pulsars 
(PSR J0106+4855) can be inferred based on the dispersion 
of the radio pulse measuring the free electron column density 
(see Section 6. 1 for details). As the remaining eight pulsars 
are radio-quiet, this method cannot be used to estimate their 
distance. Furthermore, none of the pulsars is associated with 
a known supernova remnant, preventing us from deriving dis- 
tance estimates from such source associations. 

However, it is still possible to obtain a crude estimate of the 
distance to the new pulsars, by exploiting the observed corre- 
lation between the gamma-ray luminosity L-, and the spin- 
down power E for other gamma-ray pulsars with distance 
measures (cf. Abdo et al. 2010c). Based on this correlation 
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“pseudo gamma-ray luminosities” Lp S are derived as 

Lp S - 3.2 x 10 33 (£/10 34 erg s~') 1/2 erg s“\ (12) 


where the E values are obtained from Table 2. Assuming a 
geometrical correction factor /q = 1 for the emission cone 
(Watters et al. 2009) for all gamma-ray pulsars, the rela- 
tion L~. = 4w/q Gut)d 2 is used to convert the energy flux and 
pseudo gamma-ray luminosity into a “pseudo distance” d p% , 
following Equation (2) of Saz Parkinson et al. (2010): 

03) 


1 . 6 - 


(£7 1 0 34 erg s ')^ 4 
(Gioo/lO -11 ergcnr 2 s" 1 ) 


7? k P c - 


For each of the new pulsars the resulting values forL ps and c/ ps 
are shown in Table 4. Note that these estimated gamma-ray 
luminosities and distances are subject to a number of caveats, 
detailed in Saz Parkinson et al. (2010), and could differ sig- 
nificantly from the actual values. 


5,4. Why were the new pulsars not found in previous blind 
searches ? 

To examine whether the nine pulsars found with this new 
method could be detected with previous methods, we apply 
the same search method (Atwood et al. 2006; Ziegler et al. 
2008) used to successfully discover the 26 previously-found 
blind-search LAT gamma-ray pulsars (Abdo et al. 2009a; Saz 
Parkinson et al. 2010, 2011). We select input data as done in 
previous searches; events are selected based on a fixed ROI 
and minimum-energy cut as described in Saz Parkinson et al. 
(2010). We use the same coherence window size ( T = 2 t9 s) as 
in the previous (and in the first stage of this paper’s search ) 
search. No photon weights are computed or used. No sky- 
gridding is done in the first stage of the search: only the 
2FGL-catalog sky position is used. 

The previous blind search code recovers three of the nine 
pulsars (PSRs J1620-4927, J2028+3332 and 12111+4716) 
when the 2FGL sky locations are searched. Two of the re- 
maining six pulsars are detected only if the correct known sky 
position is used (as opposed to the 2FGL-catalog position). If 
in addition the ROI and energy cuts are optimized (scanning 
different values) then three of the remaining four pulsars are 
detected. Finally, if the coherence window size is doubled 
(which dramatically raises the computational burden in a full 
blind search) then the last pulsar is found. 

As compared to the previously-published blind-pulsar 
searches of LAT data, the methods used for this work in- 
corporate several significant improvements in sensitivity and 
computational efficiency, as well as a longer data set, that ex- 
plain why these new discoveries are made. First, the use of 
efficient parameter space gridding over both sky position and 
frequency derivative allows pulsars to be found that are much 
farther from the LAT catalog position than is possible with 
previous searches. In addition, using photon weights for both 
event selection and for the search computations ensures that 
the detection significance is near optimal with only a single 
trial over event selections. In contrast, methods that use a 
"cookie cutter" event selection must either search over two 
additional parameters (minimum energy and radius for the se- 
lection), or suffer a sensitivity penalty from potentially non- 
optimal event selections. A key factor is that while weighting 
the photons provides only a modest sensitivity improvement 
over optimal cookie cutter selections, the improvement can be 
large when compared to non-optimal cuts. 


One might think that the additional systems found in this 
paper come about purely from the significantly-increased 
computer power that was available. This is not the case: the 
improved methods deserve almost all of the credit. The im- 
proved methods are so computationally efficient that had we 
searched only up to 64 Hz as was done in previous searches, 
we could have searched all 109 selected 2FGL sources over 
the initial spin-down range (|/| < 5x 10“ 13 Hz s" 1 ) and would 
have found seven of the nine new pulsars using less than 5000 
CPU hours. These are modest resources in comparison with 
those used in previous blind searches 32 . The remaining two 
pulsars would only have been found if the spin-down range 
were increased by a factor of twenty to |/[ < 10 -11 Hz s -1 , 
increasing the required CPU time to about 100000 hours. 

6. RADIO AND X-RAY COUNTERPART SEARCHES 
6. 1 . Pulsed radio emission 

These discoveries represent a substantial increase in the 
number of gamma-ray pulsars detected in blind searches of 
Fenni-LAT data. Of the 26 previously-discovered LAT blind- 
search pulsars, only 3 have been found to pulse in the radio 
band (Camilo et al. 2009; Abdo et al. 2010b; Saz Parkinson 
et al. 2010, 2011) and there are tight upper limits on the others 
(Ray et al. 201 1; Keith et al. 2011). In order to exploit these 
new discoveries in population studies of the relative beaming 
fraction and geometry of the radio and gamma-ray emission, 
it is essential to determine if they are also visible as radio pul- 
sars. 

Because the source list was chosen from Fermi- LAT pulsar- 
like unassociated sources, all of these sources have been 
previously searched for radio pulsations by the Fermi Pul- 
sar Search Consortium (Ransom et al. 2011; Hessels et al. 
2011). We have re-analyzed all of these archival observa- 
tions, with increased sensitivity because we can now do a sin- 
gle frequency trial folding the radio data with the gamma-ray 
ephemeris, and search only in a single parameter, the disper- 
sion measure (DM). This greatly reduces the number of points 
searched in parameter space (the “trials factor”) and implies 
that much smaller pulsed signal amplitudes are statistically 
significant. Where we saw an opportunity to go significantly 
deeper, we also made a number of new radio observations. 

The telescopes and observing configurations used are de- 
scribed in Table 5, and the characteristics of the individual 
radio observations are shown in Table 6. We compute the 
sensitivities using the modified radiometer equation given in 
Equation (A 1 .22) of Lorimer & Kramer (2005): 

5 -5 5T ^ I W (14) 

where 0 is the instrument-dependent factor due to digitiza- 
tion and other effects (when unknown, we assume B = 1 .25), a 
value of 5 has been assumed for the threshold signal -to-noise 
ratio for a confident pulsar detection, T sys = T iec + T s ^ y , G is 
the telescope gain, n p is the number of polarizations used (2 
in all cases), t- M is the integration time, AF is the observa- 
tion bandwidth, / is the pulsar spin frequency, and W is the 
pulse width (for uniformity, we assume W = 0.1 //). Note that 
conventionally the 3K background temperature is included in 
the receiver temperature T rec , which is measured on cold sky, 

32 The UCSC group has about 30000 CPU hours/month available: their 
previous search of 2 years of Fermi-LPtT data used 65 000 CPU hours. 
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and so T & v represents the excess temperature from the Galac- 
tic synchrotron component, which we estimate by scaling the 
408 MHz map of Haslani et al. (1982) to the observing fre- 
quency v with a spectral index a = -2.6 (defined as S „ oc v a ). 

We use a simple approximation of a telescope beam re- 
sponse to adjust the flux sensitivity in cases where the point- 
ing direction was offset from the true direction to the pulsar. 
This factor is given by q = / hwhm r / i .5 ^ w h ere q j s t h e 0 ff_ 

set from the beam center and HWHM is the beam half-width 
at half maximum. A computed flux density limit of iS Trim at 
the beam center is thus corrected via division by q for targets 
offset from the pointing direction. 

For eight of the nine pulsars, we have established that they 
are indeed radio-quiet (or extremely radio-faint), as viewed 
from Earth. 

For PSR JO 106+4855 our observations and analysis have 
revealed very faint radio pulsations. In two of our archival 45- 
minute GBT observations at 820 MHz, we detect the pulsar at 
a DM of 70.87 ± 0.2 pc cm -3 with a single-trial significance 
of 4-5 a in each observation. Although the detections are 
not individually very strong, we gain additional confidence in 
their veracity from the fact that the peak phases are consistent 
to 0.01 pulse periods when each observation is folded using 
the ephemeris determined from our gamma-ray timing. The 
summed radio profile shows a narrow peak with a duty cycle 
of ~ 2%, as shown in Figure 13. We estimate a flux density 
of ~ 20 / Jy in both observations, using the standard radiome- 
ter equation and a measurement of the off-pulse noise level. 
These radio pulsations are detected at a flux density below the 
nominal detection limit of 30 /Jy because the duty cycle is a 
factor of 5 smaller than the 10% duty cycle used in the sensi- 
tivity calculation for unknown pulse shapes. This corresponds 
to an equivalent flux density of 8 /Jy at 1400 MHz using a 
typical pulsar spectral index a = -1.7. As seen in Table 6, 
we have made two other observations of this source, one at 
350 MHz with the Green Bank Telescope and one at 1.4 GHz 
with Effelsberg. Neither of these observations detect the ra- 
dio pulsations. Accounting for this narrow duty cycle of the 
pulse, the minimum detectable flux density for the Effelsberg 
observation at 1 .4 GHz was about 1 5 //Jy, so the non-detection 
in that observation is not surprising. This constrains the spec- 
tral index to be a < -0.5. On the other hand, a spectral index 
of -1.7 would imply a flux density at 350 MHz of 85 //Jy, 
which is above the nominal sensitivity of our 350 MHz ob- 
servation. This implies that either a > -1.3 or the sensitivity 
of that observation was affected by scatter or DM broadening 
or higher than expected sky background. Using the NE2001 
model (Cordes & Lazio 2002), the measured DM corresponds 
to an estimated distance of 3.0 kpc. over a factor of 2 larger 
than the pseudo-distance. Given the radio detection, we can 
measure the so-called phase lag S = 0.073 ± 0.003 between 
the gamma-ray and radio emission. 

6.2. Continuum radio emission 

Pulsars and SNRs (Supernova Remnants) have the same 
origin, although the comparatively short lifetimes of SNRs 
means that the number of pulsar-SNR associations is quite 
small. However the rare associations are of high interest, as 
they constrain a number of pulsar and SNR parameters. A 
recent text-book example is the association of the gamma- 
ray pulsar J0007+7303 with the shell-type SNR CTA 1 (Abdo 
et al. 2008; Sun et al. 2011). 

For the nine new gamma-ray pulsars no association with 


a known SNR listed in the most recent SNR-catalog (Green 
2009) is found. However, SNRs are difficult to identify in case 
they are faint, confused or distorted by interaction with dense 
clouds. We have used large-scale and Galactic-plane radio 
continuum surveys to search for structures in the vicinity of 
the new pulsars, which may indicate an association. 

For the two pulsars PSRs JO 106+4855 and J0622+3749 lo- 
cated above 10° of Galactic latitude we used the MPIfR- 
survey sampler’ 3 to extract maps from the 408 MHz and 
1420 MHz surveys (Haslam et al. 1982; Reich 1982), which 
show faint sources and extended smooth diffuse emission, but 
no discrete object. The high-resolution interferometric NVSS 
maps (Condon et al. 1998) show numerous compact sources, 
but no extended features. 

The remaining pulsars are located within 4° of Galactic lat- 
itude covered by Galactic-plane surveys. The southern-sky 
pulsar PSR J 1620-4927 is located towards the gradient of an 
extended emission complex as seen in southern-sky surveys 
(Reich et al. 2001; Jonas et al. 1998). The Southern Galac- 
tic Plane Survey (SGPS, Haverkom et al. 2006) is insensi- 
tive to large scale emission, but shows no small scale struc- 
tures within Of 5 of the pulsar. PSR J 1746-3239 is located at 
b = -2f 2 close to a emission ridge sticking out from the plane. 
The shell-type SNR G356.6-1.5 with a size of 20' x 15' (Gray 
1994) is part of this ridge and about 0.°6 away from the pul- 
sar. PSR J 1803-2 149 is located just 2' away from the peak 
of a flat spectrum 6 Jy thermal source visible in all Galactic- 
plane surveys. Quireza et al. (2006) list details for the HII- 
region with a deconvolved Gaussian size of 2(7 and a dis- 
tance of 3.4 kpc. Also CO-emission was observed (Scoville 
et al. 1987). The pulsar is located at the periphery of this 
emission complex, but further studies are needed to investi- 
gate this region. PSR J2028+3332 is located at the southern 
boundary of the thermal Cygnus-X complex. The 1 1 cm sur- 
vey (Reich et al. 1984) shows patchy structures in the field 
and an extragalactic 1.8 Jy source 40' distant from the pul- 
sar. PSR J2030+4415 is seen towards the northern periphery 
of strong complex emission from the Cygnus-X region. Ded- 
icated studies are needed to find any emission associated with 
the pulsar. PSRs J21 11+4606 and J21 39+47 1 6 are both lo- 
cated in low emission areas, where the 21 cm and 1 1 cm Ef- 
felsberg Galactic-plane surveys show faint structures close to 
the noise level. 

6.3. X-ray 

Subsequent to the pulsar discoveries, we searched for 
archival X-ray observations covering the new pulsars’ sky 
locations. As listed in Table 7, we have found short 
Swift observations (3-10 ks exposure) for five of the pul- 
sars. In addition there is a 6 ks-long XMM-Newton obser- 
vation for PSR J 1620^4927 and a 23 ks-long Suzaku obser- 
vation for PSR J0106+4855. For the two remainins pulsars, 
PSRs J 1746-3239 and J2028+3332, 10 ks-long Swift obser- 
vations were carried out following the discoveries. 

These X-ray data were analyzed, and no X-ray counterparts 
were found for any of the new gamma-ray pulsars. 

In order to estimate upper limits on the X-ray flux 
for each new pulsar, a power-law spectrum with a pho- 
ton index of 2 and a signal-to-noise of 3 is used. The 
values for the absorbing columns are estimated analo- 
gously to Marelli et al. (2011). In the 0.3-10 keV en- 
ergy range, the derived upper limits on the X-ray flux of 

http //www.mpif r . de/ survey . html 
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the pulsars are between 1 and 3x 1CT 13 erg ernes' 1 , except 
for PSR J 1620-4927, where an upper-limit X-ray flux of 
6x 10~ 14 erg cm' 2 s _1 is obtained, and PSR 10106+4855, with 
an upper-limit X-ray flux of 2.3x 10" 14 erg cra^s" 1 . These 
results, listed in Table 7, are consistent with the gamma- 
to-X-ray flux-ratios for previously-found Fermi- LAT pulsars 
(Marelli et al. 2011). 

7. CONCLUSION 

This work reports on the discovery of nine gamma-ray pul- 
sars through the application of a new blind-search method to 
about 975 days of Fermi- LAT data. The new pulsars were 
found by searching unassociated sources with typical pulsar- 
like properties selected from the Fermi-LKT Second Source 
Catalog (Abdo et al. 201 lb). 

The sensitivity of blind searches for previously unknown 
gamma-ray pulsars is limited by finite computational re- 
sources. Thus efficient search strategies are necessary max- 
imizing the overall search sensitivity at fixed computing 
cost. We have developed a novel hierarchical search method, 
adapted from an optimal semi-coherent method (Pletsch & 
Allen 2009) together with a sliding coherence window tech- 
nique (Pletsch 2011) originally developed for detection of 
continuous gravitational-wave signals from rapidly-spinning 
isolated neutron stars. The first stage of the method is semi- 
coherent, because coherent power computed using a window 
of 6 days is incoherently combined by sliding the window 
over the entire data set. In a subsequent stage, significant 
semi-coherent candidates are automatically followed up via 
a more sensitive fully-coherent analysis. The method ex- 
tends the pioneering methods first described in Atwood et al. 
(2006). 

The new method is designed to find isolated pulsars up 
to kHz spin frequencies, by scanning a template grid in the 
four-dimensional parameter space of frequency, frequency- 
derivative and sky location. The optimal and adaptive sky 
gridding (not done in previously-published blind searches) is 
necessary, particularly at the higher search frequencies, be- 
cause the source-catalog sky-positions are not precise enough 
to retain most of the signal-to-noise ratio in year-long data 
sets. A fundamental new element of the method is the 
exploitation of a parameter-space metric (well-studied in 
the continuous gravitational-wave context (see e.g. Brady 
& Creighton 2000; Prix 2007)) to build an efficient tem- 
plate grid, as well as a metric approach to construct an op- 
timal semi-coherent combination step (Pletsch & Allen 2009; 
Pletsch 2010). A further enhancement over previous searches 
is the sub-division of the total search frequency range into 
bands via complex heterodyning. This accommodates mem- 
ory limitations, parallelizes the computational work, and per- 
mits the use of efficient sky grids adapted to the highest fre- 
quency searched in each band. A photon probability weight- 
ing scheme (Kerr 201 la) is also used for the first time in a 
published blind search, further improving the search sensitiv- 
ity. 

The nine newlv-discovered pulsars increase the previously- 
known Fermi-LKT blind-search pulsar population (Abdo et al. 
2009a; Saz Parkinson et al. 2010, 201 1) by more than one- 
third. and brings the total number to 35. The inferred pa- 
rameters of the new pulsars suggest that they belong to the 
same general population as the previously-found blind-search 
gamma-ray pulsars. Deep follow-up observations with ra- 
dio telescopes have been conducted for all of the new pul- 
sars, but significant radio pulsations have only been found for 


PSR JO 106+4855. The null results for the other eight pulsars 
indicate that they belong to the growing population of radio- 
quiet gamma-ray pulsars, which can only be detected via their 
gamma-ray pulsations. 

The computational work of the search has been done on 
the 6720-CPU-core Atlas Computing Cluster (Aulbert & 
Fehrmann 2008) at the Albert Einstein Institute in Hannover. 
Recently, in August 2011, we moved the computational bur- 
den of the search onto the volunteer distributed computing 
system Einstein@Home i4 (Abbott et al. 2009a,b; Knispel 
et al. 2010). This will provide significantly more computing 
power, and will allow a complete search of the parameter- 
space up to kHz pulsar spin frequencies. We also hope that in 
the future an improved version of these methods can be used 
to carry out blind searches for gamma-ray pulsars in binary 
systems. 

The combination of improved search techniques and much 
more powerful computational resources leave us optimistic 
that we can find still more gamma-ray pulsars in the Fenni- 
LAT data. These advances should also greatly increase the 
chance of finding the first radio-quiet gamma-ray MSP with 
the Fermi-LKT. We hope that further discoveries, and further 
study of these systems, will eventually provide important ad- 
vances in our understanding of pulsars, and of their emission 
mechanisms and geometry. 
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TABLE 1 

Names and Sky Locations of the Discovered Gamma-Ray Pulsars 


Pulsar Name 

Source Association 

R.A. a 

(hh:mm:ss.s) 

Decl.® 

(dd:mm:ss.s) 

/ b 

(deg) 

b h 

(deg) 

J01 06+4855 

2FGL J01 06.5+4854 
1 FGL J0106. 7+4853 

01:06:25.06(1) 

+48:55:51.8(2) 

125.5 

-13.9 

J0622+3749 

2FGLJ062 1.9+3750 
1 FGL J0622. 2+3751 

06:22:10.51(2) 

+37:49:13.6(9) 

175.9 

11.0 

J 1620-4927 

2FGLJ 1620.8-4928 
1 FGL J 1 620.8-4928c 
1AGLJ 1624-4946 

16:20:41.52(1) 

-49:27:37.1(3) 

333.9 

0.4 

J 1746-3239 

2FGLJ 1746.5-3238 
1 FGL J 1746.7-3233 

17:46:54.947(8) 

-32:39:55.8(7) 

357.0 

-2.2 

J 1803-2 149 

2FGLJ 1803.3-2 148 
1FGLJ1803.1-2147C 
1AGLJ 1805-2 143 

18:03:09.632(9) 

-21:49:13(4) 

8.1 

0.2 

J2028+3332 

2FGLJ2028. 3+3332 
3EG J2027+3429 

20:28:19.860(5) 

+33:32:04.36(7) 

73.4 

-3.0 

J2030+4415 

2FGL J2030.7+4417 
1FGLJ2030.9+441 1 

20:30:51.35(4) 

+44:15:38.1(4) 

82.3 

2.9 

J21 11+4606 

2FGLJ21 11.3+4605 
1 FGL J2 11 1.3+4607 
OFGL J2 110.8+4608 

21:11:24.13(3) 

+46:06:31.3(3) 

88.3 

-1.4 

J2139+4716 

2FGLJ2 139.8+47 14 
1 FGL J 2 1 39.9+47 1 5 

21:39:55.95(9) 

+47:16:13(1) 

92.6 

-4.0 


Note. — A list of the nine new pulsars reported in this work, showing their sky locations 
and associations with cataloged gamma-ray sources. The associations listed include sources from 
the Fermi - LAT Second Source Catalog (2FGL. Abdo et al. 201 lb), the Fermi - LAT First Source 
Catalog (1FGL, Abdo et al. 2010a). the Fermi-LAJ Bright Source List (OFGL. Abdo et al. 
2009b), the AGILE Catalog (1AGL, Pittori et al. 2009). and the Third EGRET Catalog (3EG, 
Hartman etal. 1999). 

a Right ascension (J2000.0) and declination (J2000.0) obtained from the timing model, where the 
numbers in parentheses are the statistical 1 a errors in the last digits. 
b Galactic longitude (/) and latitude (70. rounded to the nearest tenth of a degree. 


TABLE 2 

Measured and Derived Parameters of the Discovered Gamma-Ray Pulsars 


Pulsar Name 

/ 

(Hz) 

/ 

(-10' 13 Hz s' 1 ) 

Weighted 
FI- test 

T 

(kyr) 

E 

(10 34 erg s' 1 ) 

As 

(10 ,2 G) 

Slc 

(kG) 

JO 106+4855 

12.02540173638(8) 

0.61881(7) 

843.1 

3081.1 

2.9 

0.2 

3.0 

J0622+3749 

3.00112633651(5) 

2.28985(4) 

288.8 

207.8 

2.7 

2.9 

0.7 

J 1620-4927 

5.81616320951(5) 

3.54782(4) 

566.4 

259.9 

8.1 

1.4 

2.4 

J 1746-3239 

5.01149235750(3) 

1.64778(3) 

249.8 

482.2 

3.3 

1.2 

1.3 

J 1803-2 149 

9.4044983174(2) 

17.25894(6) 

451.9 

86.4 

64.1 

1.5 

11,0 

J2028+3332 

5.65907208453(2) 

1.55563(2) 

1108.3 

576.8 

3.5 

0.9 

1.5 

J2030+4415 

4.4039248637(5) 

1.2576(2) 

584.8 

555.2 

2 2 

1.2 

1.0 

J21 11+4606 

6.3359340865(4) 

57.4218(3) 

554.3 

17.5 

143.6 

4.8 

11.1 

J2 139+47 16 

3.5354509962(2) 

0.2232(2) 

351.1 

2511.5 

0.3 

0.7 

0.3 


Note. — The reference epoch for all measured rotational parameters is MJD 55225 and the time range for ail 
timing models is MJD 54682 - 55719. The derived quantities in columns 5-8 are based on the / and / values 
obtained from the timing model and are rounded to the nearest significant digit. To model the timing noise present 
in PSR J 1803-2149. a second frequency derivative is necessary: / =7,3(8)x 1CT 24 Hz s~ 2 . To model the timing noise 
present in PSR .12030+4415. higher frequency derivatives up to third order are necessary: / = ~1.5(3)x IO~ 23 Hz s' 2 
and f = -6(2)x 10~ 31 Hz s' 3 . To model the riming noise present in PSR J21 1 1+4606, higher frequency derivatives 
up io fourth order are necessary: /=2.30<5)x 10~ 22 Hz s~ 2 , / =-7.9(2)x 10~ 30 Hz s~ 3 and / =3.2(4) x 10~ 3, Hzs -4 . 
The numbers in parentheses are the statistical lcr errors in the last digits. 
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JO 106+4855 



Pulse phase Pulse phase 


Fig. 4. — Phase-time diagram and pulse profile for PSR J0106+4855. The 
left panel shows the pulse phase at the arrival time of each photon, where 
the gray-scale intensity represents the photon probability weight. The upper 
right plot shows the summed probability weights: the integrated pulse profile 
using a resolution of 32 bins per rotation. The error bars represent the hr 
statistical uncertainties. The four plots below resolve the integrated pulse 
profile according to separate energy ranges. For clarity, the horizontal axis 
shows two pulsar rotations in each diagram. In obtaining these plots the 8000 
events with the highest probability weights have been used (as in the blind 
search). 



Pulse phase Pulse phase 


Fig. 6. — Phase-time diagram and pulse profile for PSR J 1620-4927. The 
plots have identical form as those shown in Figure 4. Note that for this par- 
ticular pulsar, based on the photon probability weights, no selected events 
(among the 8000) have energies in the 0. 1-0.3 GeV range. 


.10622+1749 



Fig. 5. Phase-time diagram and pulse profile for PSR J0622+3749. The 

piots have identical form as those shown in Figure 4. Fig. 7. — Phase-time diagram and pulse profile for PSR J 1746-3239. The 


plots have identical form as those shown in Figure 4. 
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Fig. 8. — Phase-time diagram and pulse profile for PSR J 1803-2149. The FlG. 10. — Phase-time diagram and pulse profile for PSR J2030+4415. The 

plots have identical form as those shown in Figure 4. plots have identical form as those shown in Figure 4. 




Pulse phase 


J2 1 1 1+4606 




Fig. 9.— Phase-time diagram and pulse profile for PSR J2028+3332. The FIG. 1 1 .— Phase-time diaaram and pulse profile for PSR J21 1 1+4606. The 
plots have identical form as those shown in Figure 4. plots have identical form as those shown in Figure 4. 
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J2i3<M47!6 




20 



Pul se phase 


Pulse phase 


Fig. 12. — Phase-time diagram and pulse profile for PSR J2139+4716. The 
plots have identical form as those shown in Figure 4. 
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Pulse phase 

FIG. 13. — Phase-aligned gamma-ray (top) and radio (bottom) pulse pro- 
files for PSR J0106+4855. 
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TABLE 3 

Pulse-Profile Parameters of the Discovered Gamma-Ray Pulsars 


Pulsar Name Peak Multiplicity FWHM[ 

FWHM; 


A 


J0106+4855 

2 

0.02 

± 

0.01 

0.05 

rfc 

0.01 

0.50 

± 

0.01 

J0622+3749 

2 

0.21 

± 

0.04 

0.03 

± 

0.03 

0.47 

± 

0.03 

J 1620-4927 

2 

0.13 

± 

0.04 

0.15 

± 

0.05 

0.21 

± 

0.03 

J 1746-3239 

2 

0.32 

± 

0.07 

0.11 

± 

0.11 

0.63 

± 

0.06 

J 1803-2 149 

2 

0.04 

± 

0.02 

0.05 

± 

0.03 

0.40 

± 

0.02 

J2028+3332 

2 

0.10 

± 

0.02 

0.04 

± 

0.01 

0.38 

± 

0.02 

J2030+4415 

2 

0.09 

± 

0.03 

0.04 

± 

0.01 

0.49 

± 

0.02 

J21 11+4606 

2 

0.08 

± 

0.02 

0.07 

± 

0.02 

0.31 

± 

0.02 

12139+4716 

1 

0.13 

4- 

0.03 








Note. — For each of the nine pulsars, we give the parameters describing 
the shape of the pulse profile, including the peak multiplicity, the Full- Widths at 
Half Maxima (FWHM) of the peaks, and the separation A between the gamma-ray 
peaks for pulsars with more than one peak. 


TABLE 4 

Spectral Parameters of the Discovered Gamma-Ray Pulsars 


Pulsar Name 

r 

E c 

W 

G»X) b 

^-ps ^ps 



(GeV) 

(10' 8 photons cm' 2 s' 1 ) 

( 1 0 -1 1 erg cm -2 s' 1 ) 

(10 33 ergs _l ) (kpc) 


JO 106+4855 

1.47 

± 

0.23 

± 

0.12 

3.31 

± 

0.92 

± 

0.08 

2.56 

± 

0.77 

±0.32 

2.40 ± 

0.31 

±0.04 

5.5 

1.4 C 

J0622+3749 

0,59 

± 

0.34 

± 

0.09 

0.60 

± 

0.13 

± 

0.04 

2.21 

± 

0.35 

±0.08 

1.69 ± 

0.15 

±0.05 

5.3 

1.6 

J 1620-4927 

1.01 

± 

0.18 

± 

0.05 

2.44 

± 

0.42 

± 

0.28 

9.61 

± 

1.68 

±0.94 

13.5 ± 1 

,0± 

1.7 

9.1 

0.7 

J 1746-3239 

1.33 

± 

0.08 

± 

0.35 

1.65 

± 

0.12 

± 

0.52 

9.97 

4- 

0.94 

± 1.85 

7.86 ± 

0.41 

hh 

o 

-4 

5.8 

0.8 

J 1803-2 149 

1.96 

± 

0.11 

± 

0.20 

5.73 

± 

1.72 

± 

2.07 

20.7 ± 3 

1.1 ± 

0.5 

13.1 ± 1 

.1 ± 

2.1 

25.6 

1.3 

J2028+3332 

0.86 

± 

0.21 

± 

0.07 

1.53 

± 

0.24 

± 

0.08 

5.12 

± 

0.87 

±0.44 

6.09 ± 

0.41 

±0.13 

6.0 

0.9 

J2030+4415 

1.89 

± 

0.14 

± 

0.22 

2.16 

± 

0.65 

dt 

0.67 

13.3 ± 1 

.4± 

0.2 

7.06 ± 

0.48 

±0.66 

4.8 

0.7 

J21 1 1+4606 

1.63 

± 

0.14 

± 

0.05 

5.43 

± 

1.80 

± 

1.56 

4.39 

± 

0.69 

± 0.02 

4.13 ± 

0.34 

±0.30 

38.4 

2.7 

J2 139+47 16 

0.80 

± 

0.27 

± 

0.02 

1.02 

± 

0.21 

± 

0.07 

2.65 

4- 

0.44 

±0.19 

2.51 ± 

0.21 

±0.01 

1.8 

0.8 


Note. — This table describes the spectral properties of each of the nine pulsars, modeling each spectrum as an exponentially cut-off 
power-law with photon indices T and cutoff energies E c . The spectral parameters listed here for each pulsar are obtained from maxi- 
mum likelihood fits. The first quoted uncertainties are statistical, while the second are systematic and correspond to the differences in the 
best-fit parameters observed when doing the spectral analyses with the P6„V3 iRFs and associated diffuse emission models (namely, the 
g ll_iem_v02 map cube and isotropic _iem_v02 template). For each object, the pseudo gamma-ray luminosity L ps and the pseudo distance d ps 
are inferred from the apparent spin-down power E and the energy flux Gim above 100 MeV. Note that these estimated gamma-ray lumi- 
nosities and distances are subject to a number of caveats, detailed in Saz Parkinson et al. (2010), and could differ significantly from the 
actual values. 

a Photon flux measured above 100 MeV. 
b Energy flux measured above 100 MeV. 

c The actual distance is 3.0 kpc, as inferred from the dispersion of the radio pulse measuring the free electron column density; see Section 6. 1 . 


TABLE 5 

Definition of Radio Observing Codes 


Obs Code 

Telescope 

Gain 

(K/Jy) 

Frequency 

(MHz) 

Bandwidth A F 
(MHz) 

/3 a 

"p 

HWHM 

(arcmin) 

Tree 

(K) 

GBT-350 

GBT 

2.0 

350 

100 

1.05 

2 

18.5 

46 

GRT-820 

GBT 

2.0 

820 

200 

1.05 

2 

7.9 

29 

GBT-S 

GBT 

1.9 

2000 

~0(P 

1.05 

2 

3.1 

22 

Eff-Ll 

Effelsberg 

1.5 

1400 

250 

1.05 

2 

9.1 

22 

E1T-L2 

Effelsberg 

1.5 

1400 

140 

1.05 

2 

9.1 

22 

Jodrell 

Lovell 

0.9 

1520 

200 

1.05 

2 

6.0 

24 

AO-327 

Arecibo 

11 

327 

25 

1.12 

7 

6.3 

116 

AO-Lwide 

Arecibo 

10 

1510 

300 

1.12 

2 

1.5 

27 

Parkes-BPSR 

Parkes 

0.735 

1352 

340 

1,05 

9 

7.0 

25 


Note. — The sky locations of all nine pulsars have been searched for pulsating radio 
emissions. This Table gives the radio telescope and back-end parameters used in those 
observations, which are described in Table 6. 


a Instrument-dependent sensitivity degradation factor. 

b The instrument records 800 MHz of bandwidth, but to account for a notch filter for RFI 
and the lower sensitivity near the band edges, we use an effective bandwidth of 700 MHz 
for the sensitivity calculations. 
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TABLE 6 

Radio Search Observations of the New Gamma-Ray Pulsars 


Target 

Obs Code 

Dale 

Ant 

(min) 

R.A. a 

(J2000) 

Deck 3 

(J2000) 

Offset 

(arcmin) 

T s ky 

(K) 

^min 

y ) 

JO 106+4855 

GBT-350 

2009-10-25 

32 

01:06:37.7 

48:54:11 

2.7 

49.3 

136 


GBT-820 

2010-11-17 

45 

01:06:35.5 

48:55:30 

1.8 

5.4 

30 b 


GBT-820 

2010-12-17 

45 

01:06:35.5 

48:55:30 

1.8 

5.4 

30 b 


Eff-L2 

2011-06-01 

45 

01:06:25.1 

48:55:52 

0.0 

1.3 

31 

J0622+3749 

GBT-350 

2009-10-27 

32 

06:22:05.5 

37:51:07 

2.2 

46.2 

131 


Eff-Ll 

2010-05-14 

32 

06:22:14.7 

37:51:49 

2.8 

1.3 

30 


Eff-Ll 

2010-02-06 

10 

06:22:15.0 

37:51:48 

2.8 

1.3 

53 


Eff-Ll 

2010-07-10 

52 

06:22:13.0 

37:50:36 

1.5 

1.3 

22 


Eff-Ll 

2010-07-10 

55 

06:22:13.0 

37:50:36 

1.5 

1.3 

22 


GBT-820 

2010-12-12 

45 

06:21:59.0 

37:51:36 

3.3 

5.0 

32 


GBT-820 

2010-12-17 

45 

06:21:59.0 

37:51:36 

3.3 

5.0 

32 

J 1620-4927 

Parkes-BPSR 

2009-08-03 

270 

16:21:05.5 

-49:30:32 

4.9 

16.9 

42 


Parkes-BPSR 

2010-11-18 

144 

16:20:43.5 

-49:28:24 

0.9 

16.9 

42 


Parkes-BPSR 

2011-05-10 

72 

16:20:41.3 

-49:27:36 

0.0 

16.9 

58 

J 1746-3239 

GBT-S 

2009-12-23 

60 

17:46:47.9 

-32:36:22 

3.8 

5.1 

30 


GBT-820 

2010-11-14 

45 

17:46:41.0 

-32:36:18 

4.6 

51.6 

85 


Parkes-BPSR 

2011-05-10 

72 

17:46:54.9 

-32:39:55 

0.0 

14.1 

54 

J 1803-2 149 

Eff-Ll 

2010-02-13 

25 

18:03:12.0 

-21:47:27 

1.6 

17.8 

55 


Eff-Ll 

2010-05-22 

32 

18:03:11.7 

-21:47:28 

1.5 

17.8 

48 


GBT-S 

2010-09-04 

65 

18:03:11.7 

-21:47:28 

1.5 

7.1 

13 

J2028+3332 

GBT-820 

2009-08-13 

60 

20:27:48.0 

33:32:24 

6.6 

10.3 

46 


GBT-S 

2010-09-20 

30 

20:28:18.0 

33:33:23 

1.4 

1.0 

15 


GBT-820 

2010-11-22 

45 

20:28:19.0 

33:32:53 

0.8 

10.3 

33 


GBT-820 

2010-12-17 

45 

20:28:19.0 

33:32:53 

0.8 

10.3 

33 


AO-Lwide 

2011-05-21 

45 

20:28:19.9 

33:32:06 

0.0 

2.1 

4 


AO-327 

2011-05-30 

25 

20:28:19.9 

33:32:06 

0.0 

112.9 

142 

J2030+4415 

Eff-Ll 

2010-02-07 

10 

20:30:55.0 

44:11:52 

3.8 

4.0 

62 


Eff-Ll 

2010-05-15 

32 

20:30:55.3 

44:11:53 

3.8 

4.0 

35 


Eff-Ll 

2010-07-10 

60 

20:30:59,2 

44:15:33 

1.4 

4.0 

23 


Eff-Ll 

2010-07-30 

60 

20:30:59.2 

44:15:33 

1.4 

4.0 

23 


GBT-820 

2010-11-22 

45 

20:30:54.7 

44:16:08 

0.8 

16.0 

38 


GBT-820 

2011-05-28 

183 

20:30:51.3 

44:15:38 

0.0 

16.0 

19 


Eff-L2 

2011-06-01 

45 

20:30:51.5 

44:15:37 

0.0 

4.0 

35 

J21 11+4606 

GBT-820 

2009-09-19 

60 

21:11:22.8 

46:05:53 

0.3 

16.0 

33 


Jodrell 

201 l-06-22 c 

14x60 

21:11:24.0 

46:06:29 

0.0 

9.0 

14 

J2 139+47 16 

Eff-Ll 

2010-07-10 

60 

21:39:52.3 

47:13:43 

2.6 

2.0 

22 


GBT-350 

2009-10-25 

32 

21:39:53.2 

47:15:22 

1.0 

74.9 

171 


GBT-820 

2010-12-11 

45 

21:39:53.5 

47:13:30 

2.7 

8.2 

34 


GBT-820 

2010-12-18 

45 

21:39:53.5 

47:13:30 

2.7 

8.2 

34 


Eff-Ll 

2010-05-15 

32 

21:39:56.9 

47:15:28 

0.8 

2.0 

29 


Note. — The sky locations of all nine pulsars have been searched for radio pulsations. Only for 
PSR JO 106+4855 radio pulsations are detected. The minimum detectable flux density <Smin for each 
observation is computed at the observing frequency using Equation (14) and the parameters in Table 
5. as described in the text. 

a Telescope pointing direction (not necessarily source position) 

b For these two observations, radio pulsations were detected at a flux density of 20 /LtJy (see text). 
c Observed 14 times for 1 hour each between this date and 201 1-07-1 1 . 
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TABLE 7 

X-ray Coverage of the Discovered Gamma-Ray Pulsars 


Pulsar Name 

Instrument 

Exposure Time 
(ks) 

Absorbing column 11 
(10 21 cm -2 ) 

Flux Upper Limit b 
(10~ 13 erg cm' 2 s" 1 ) 

L-,/Ty c 

JO 106+4855 

Suzaku 

23.0 

1.0 

0.843 

>285 

J0622+3749 

Swift XRT 

4.4 

1.0 

2.58 

>65 

J 1620-4927 

XMM-Newton 

6.0 

4.0 

0.674 

>1988 

J 1746-3239 

Swift XRT 

8.7 

1.0 

1.74 

>451 

J 1803-2 149 

Swift XRT 

7.7 

5.0 

3.26 

>34 

J2028+3332 

Swift XRT 

10.3 

1.0 

1.57 

>387 

J2030+4415 

Swift XRT 

10.2 

4.0 

2.53 

>279 

J21 11+4606 

Swift XRT 

10.1 

3.0 

2.25 

>183 

J2139+4716 

Swift XRT 

3.2 

1.0 

3.20 

>78 


Note. — The sky locations of all nine pulsars have been searched for (non-pulsating) X-rays, 
using both archival and new data. No X-ray sources were found at the new pulsar locations, so 
flux upper limits and lower limits on the gamma-ray to X-ray luminosity ratio are reported. 
a Estimated analogously to Marelli et al. (201 1). 

b Upper limit on the unabsorbed flux in the 0.3- 10 keV energy range, using an absorbed power- 
law model with a photon index of 2 and a signal-to-noise of 3. 
c Lower limit on the L-,/Lx flux ratio. 
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